Quantum Simulation of Tunneling in Small Systems 
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A number of quantum algorithms have been performed on small quantum computers; these include 
Shor's prime factorization algorithm, error correction, Grover's search algorithm and a number of 
analog and digital quantum simulations. Because of the number of gates and qubits necessary, 
however, digital quantum particle simulations remain untested. A contributing factor to the system 
size required is the number of ancillary qubits needed to implement matrix exponentials of the 
potential operator. Here, we show that a set of tunneling problems may be investigated with no 
ancillary qubits and a cost of one single-qubit operator per time step for the potential evolution. 
We show that physically interesting simulations of tunneling using 2 qubits (i.e. on 4 lattice point 
grids) may be performed with 40 single and two-qubit gates. Approximately 70 to 140 gates are 
needed to see interesting tunneling dynamics in three-qubit (8 lattice point) simulations. 

PACS numbers: 



Quantum simulations on quantum computers are one of 
a set of algorithms that give exponential improvement in 
computational resources relative to the best classical al- 
gorithm. Small quantum simulations have already been 
realized on NMR JM, atomic @, [nj, ion trap |ll|^ 
and photonic lg, [l7[ quantum computers in essentially 
two forms: 1) analog simulations in which a quantum 
Hamiltonian (typically many-body or multiple spin) is 
mapped either directly or via a suitable pulse-sequence to 
a computational Hamiltonian, and 2) digital simulations 
in which a quantum system's Hamiltonian is split into 
free and interacting operators, then, using Trotter's for- 
mula, is simulated on a quantum computer. Digital quan- 
tum particle simulations [3] such as those proposed for 
the simulation of chemical dynamics [l9| have remained 
untested due to the large number of gates and/or ancil- 
lary qubits needed to compute the kinetic and potential 
operators IjJ 20j]. 

Digital Quantum Particle Simulations The standard 
digital quantum simulation alg orithm for a particle on 
a one-dimensional grid 



ficiently in 



log 2 A 



18l . |20j encodes the position ef- 
qubits, where N is the size of 



the lattice of discretized particle locations, x% = kAx, 
k = 0, . . . , N — 1. The method uses a split operator ap- 
proach to integrate a Schrodinger equation with a time- 
independent Hamiltonian that is first-order accurate in 
the time step, At [Hj]: 
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Higher order methods that give more accurate time in- 
tegration have been developed 21- 2^, but methods of 
order higher than two require more gates per time step. 
Therefore we will only consider first-order methods here. 
States in the qubit Hilbert space 



where \ipi) € {|0), |1)}, represent particle location in the 



binary representation, \x) = 



j'eo, 



The 



matrix exponential for the kinetic operator is calculated 
using a quantum Fourier transform (QFT) 

e -iKt = Fe -iTtp\ j (2) 

where is a discrete Fourier transform operator and 
T is diagonal with entries proportional to — q 2 /2m, and 
q denotes the Fourier mode wavcnumbcr. The resulting 
periodic, shift invariant unitary transform, e~ lKAt , gives 
an approximation to the time evolution on the lattice due 
to the kinetic energy operator that is accurate to N'th 
order in space. This leads to the digital quantum particle 
simulation algorithm: 
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The QFT takes of order n 2 gates to calculate [2J| and 
general algorithms implementing the diagonal T and V 
operators require ancillary qubits (25j . although it has 
been shown that the quadratic kinetic energy operator 
may be computed with n 2 two-qubit gates with no ancil- 
lary qubits [20l | . Thus, the qubit count for this algorithm 
is dominated by the calculation of the diagonal potential 
operator. 

Square-well Potentials A special case of considerable 
interest is that of the square well potential. Square-well 
potentials are commonly used in the study of a number 
of quantum tunneling phenomena. A set of square-well 
potentials may be implemented with a sole single-qubit 
operator and no ancillary qubits. To see this, consider 
the single-qubit Z-rotation on the highest order qubit 
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where v is a parameter, a superscript indicates the qubit 
to which the operator is applied and a z is the Pauli z- 
matrix 
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FIG. 1: Particle probability distributions as a function of time 
for the first four steps of a two qubit simulation for v = (free 
particle) and v = 10 (particle in double well). The double- 
well potential and gray scale used to plot the probabilities are 
depicted to the right. 



The operator e~ lVAt , when acting on a lattice state, im- 
plements a square- well potential by rotating qubit states 
with |0) (|1) resp.) highest order qubit with positive (neg- 
ative resp.) phase velocity v. The single-qubit operator 
acting on the next highest order qubit 
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implements a double-square-well potential, and so on, 
with the last potential in this series implementing a 
comb-like potential. 

A Two Qubit Tunneling Simulation Let us consider the 
smallest possible tunneling simulation. An n = 2 qubit 
simulation of N = 4 lattice points may be performed with 
a double well potential, exp(— iVAt) = exp(— i ya^ At). 
The QFT may be computed with the operators [24j 



F' — HqQqiHi 
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where Hi is a Hadamard operator on the i'th qubit 
and the controlled-phase gate fioi = diag(l, 1, l,u>) with 
u> = exp(27ri/4). Note that results in a bit-swapped 
Fourier transform (in our notation, F is the inverse 
Fourier transform matrix). 

The kinetic operator is then K = FDF\ with D = 
exp(-i (-2tt/4) 2 diag(0,4, 1, l)At) (Note that this oper- 
ator is also bit-swapped and we have taken m = 1/2). 
This operation may be achieved up to an overall phase 
with the operator 



D = QoiZ^o 



where the single-qubit operators 
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and $oi is a controlled-phase operator on qubits and 
1, 



01 



_ p -il c 2 diag(l,l,l,-l)oiAt 



(9) 



where 7 = (— 2tt / 4) 2 / The coefficients (cq = — 1, 
ci = —4 and C2 = 4) in the unitary operators Zq, Z\ and 
$01 (resp.) were obtained by noting that the vectors 
(1,1,-1,-1), (1,-1,1,-1) and (1,1,1,-1) that form 
their diagonal elements are a basis for zero mean vec- 
tors (i.e. neglecting the constant phase proportional to 
(1,1,1,1)) for R 4 . 
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One Time Step of a Two-qubit Digital Quantum Single-particle Simulation Circuit 

I 



This method of forming a diagonal Hamiltonian using a 
basis of diagonal operators is not, in general, efficient and 
scales as 2™ — 1. However, for small systems of qubits it 
takes fewer gates. For instance, the Benenti-Strini (BS) 
algorithm [20| for forming a quadratic function on the 
diagonal requires n 2 two-qubit operators. Note, how- 
ever, that for n = 2,3 and 4, we would need 4,9 and 
16 (resp.) operators for BS, but at most 3, 7 and 15 
(resp.) are required using a basis along the diagonal. 
Furthermore, fewer multi-qubit operations are necessary 
with this method, although a three-qubit operator may 
be necessary for three-qubit simulations and three- and 
four-qubit operators for four-qubit simulations. Here, we 
use the method to construct the diagonal operator D for 



the kinetic energy operator, however, arbitrary potentials 
could also be constructed this way. Because our goal is 
to find interesting simulations with few gates, we do not 
pursue more complex potentials here. 

In total, each time step in this tunneling simulation 
requires 10 operations, 7 single qubit operations and 3 
two-qubit operations. The circuit for a single time step 
is shown above. This circuit implements a double-well 
potential with the gate Po = exp(— iva^At) acting on 
the lowest-order qubit. 

In Fig. 1, we plot lattice occupation probabilities 
from two simulations with a double-well potential with 
At = 1/10: one a free-particle simulation with v = 
and the other a tunneling simulation with v = 10. The 
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initial state was \ipi n it) = 1 01), corresponding to a par- 
ticle in one of the wells. The free-particle probability 
distribution spreads across all lattice points as it evolves, 
whereas the particle tunnels from the well at lattice point 
1 (| 01)) to the well at lattice point 3 (|11)) in the tun- 
neling simulation. These results show that differences in 
the evolution of the probability distribution are evident 
within 4 time steps. Thus, such a simulation may be 
implemented on a quantum computer with 4 x 10 = 40 
gates (28 single-qubit and 12 two-qubit). 

Multi-Qubit Quantum Tunneling Simulations Larger 
simulations require more gates. For instance, a three- 
qubit simulation requires 6 gates per QFT (3 single-qubit 
and 3 two-qubit), 6 gates for the diagonal kinetic en- 
ergy operator (3 single-qubit and 3 two-qubit), and one 
single-qubit gate for the potential as shown in the circuit 
diagram below. 
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where Pi = exp(-ivalAt) is shown, representatively, for 
a double-well potential, but other square-well potentials 
could be generated by acting on different qubits, 
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Here, the unitary operators are the two-qubit controlled- 
phase operators Clij = Qij(uj), where co = cxp(27rz/8), 
the single-qubit operators 



Z 
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exp (— i'jcoa^At) 
exp (—i / ycia z At) 
exp (—ijc 2 o-lAt) 



and the two-qubit controlled-phase operators 

$ 01 = exp(-i7c 3 diag(l,l,l,-l) iAt) 
$02 = exp(-i7c 4 diag(l, 1, 1, -l) 02 Ai) 
$12 = exp (-Z7c 5 diag(l, 1, 1, -l) 12 At) , 



where 



and 
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c = -1.42 
ci = —5.66 
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FIG. 2: Particle probability distribution as a function of time 
over ten time steps in a three-qubit double-well simulation. 
The potential schematic is shown to the right, as in Fig. 1. 
Gray scales are as in Fig. 1. The double- well potential has 
two lattice points per well and, for the given initial state, an 
oscillation is induced in one of the wells then tunnels to the 
other well. 



C2 
C3 
C 4 
C5 



-22.63 
22.63 
11.31 
-5.66 . 



Note that, in principle, a seventh three-qubit opera- 
tor would also be necessary proportional to the vec- 
tor (1, 1, 1, 1, 1, 1, 1, —1) along the diagonal, but for the 
(bit-swapped) diagonal of the 8 lattice point simulation, 
(0,16,4,4,1,9,9,1), its coefficient is identically zero. 
This circuit requires a total of 19 gates (10 single-qubit 
and 9 two-qubit). 

In Fig. 2, we show results from a three-qubit simu- 
lation with a double-well potential, where each well is 
resolved with two lattice points. The time step At — 1/5 
and v = 5. The initial state was \4>init) = |H0). Because 
the initial state only occupies half of one well, oscillatory 
dynamics are visible within the well. After a few time 
steps, the oscillatory state tunnels between wells. Oscil- 
latory dynamics are evident within 4 time steps, but 5 
time steps are required before the oscillatory state tun- 
nels appreciably to the second well and 7 time steps are 
needed to see oscillation of the tunneled state. Thus, be- 
tween 4 x 19 = 76 or 7 x 19 = 133 gates would be required 
in order to see interesting tunneling effects in such a sim- 
ulation. Four-qubit simulations can be envisaged using 
circuits based on similar methodology. 

Conclusion Only n = log 2 N qubits are required for 
an N lattice-point particle simulation, therefore this al- 
gorithm is efficient in the number of qubit resources re- 
quired. With a very few qubits, interesting tunneling 
dynamics may be simulated with a gate count that is 
within reach of current quantum architectures. 
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